From QCD to a dynamical quark model: construction and some meson spectroscopy 
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We introduce an effective quark model that is in principle dynamically derivable from the QCD action. An 
important feature is the incorporation of spontaneous chiral symmetry breaking in a renormalizable fashion. 
The quark propagator in the condensed vacuum exhibits complex conjugate poles, indicative of an unphysical 
spectral form, i.e. confined quarks. Moreover, the ensuing mass function can be fitted well to existing lattice 
data. To validate the physical nature of the new model, we identify not only a massless pseudoscalar (i.e. a pion) 
in the chiral limit, but we also present reasonable estimates for the p meson mass and decay constant, employing 
a contact point interaction and a large N argument to simplify the diagrammatic spectral analysis. We stress that 
we do not use any experimental input to obtain our numbers, but only rely on our model and lattice quark data. 
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Next to confinement — the absence of color charge par- 
ticles as in the observable QCD spectrum, the other crucial 
nonperturbative ingredient of QCD is the dynamical breaking 
of the (almost) chiral symmetry (D^SB). The latter breaking 
is what ensures most of the mass of the mesons/baryons and 
what explains the mass gap between the light pseudoscalar 
mesons and the rest of the spectrum due to the Goldstone bo- 
son nature of the former. Albeit confinement and D^SB are 
well appreciated, their treatment, as deeply nonperturbative 
phenomena, inspires up to today active research. Here, we 
present a novel quark model, directly grounded in QCD. We 
have a few prerequisites to be met, i.e. the model should be (i) 
able to capture the correct quark dynamics, in particular to de- 
scribe properly the nonperturbative lattice quark propagator; 
(ii) renormalizable, as we require to describe quark dynamics 
over the whole momentum range, including its (well-known) 
UV behaviour; (iii) practical to compute with; (iv) display- 
ing the correct chiral behaviour such as massless pions in the 
chiral limit; (v) allowing to construct the rest of the meson 
spectrum. 

We initiate from the QCD action in Ad Euclidean space 
[22], but we add an extra piece to it: 

S = Jd 4 x ^ l -F% + y]py -l&% - \d 2 X - r|3 2 9 + 93 2 ri ) 

The additional (fermion) fields are perturbatively trivial, as the 
underlined (quadratic) piece constitutes a unity. Moreover, it 
is BRST exact [1]. The previous action is in fact equivalent 
to the QCD one, in particular should its symmetry content 
be. We treat the new fermion fields as chiral singlets, as their 
quadratic form is not of the usual kind [23]. We nevertheless 
recover the QCD chiral symmetry with the transformations 

8 5 \|/ = ;y 5 \|/, 8 5 ? = /Ws, 8 5 (rest)=0. (1) 

We did not explicitly write the gauge fixing part. As our ac- 
tion is supposed to properly describe nonperturbative gluon 
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and quark dynamics, a suitable nonperturbative gauge fixing 
as the Gribov-Zwanziger scheme in the Landau gauge can be 
selected [2]. It is now useful to record the following mass 
dimensions dim[\|/,\j7] =3/2 and dim[A,,c^,r|,8, A,,^,ff,8] = 1. 
Consider next the local composite operators (LCO) 0\ = 
vj/£ - A,\|/ — vj/A, and O2 = M; + %X + rf0 - 8r|, with dim[0i] = 
5/2 and dim[02] = 2. The mixed fermion condensate (0\) 
serves as an order parameter for chiral symmetry breaking, in- 
deed Oi = -8571 with 71 = -i (fy{ 5 \\f + \j/y 5 £ - Xj 5 \\t - \jry 5 A,) . 

We shall later on show that this 7C corresponds to the pion 
field, viz. a massless pseudoscalar in the here considered chi- 
ral limit. Let us now construct the underlying effective action 
so that the nonperturbative dynamics related to (0\) can be 
taken into due account. We introduce 2 scalar sources, J and 
j, coupled to 0\ and O2, with dim[/] = 3/2, dim[/'] = 2, 

S^S + Jd 4 x(jO l +j0 2 -^g 2 ) J -pJ (2) 

as acting with a derivative w.r.t the sources allow to define 
the operators at the quantum level. The parameter £ must be 
introduced to allow for a homogeneous linear renormaliza- 
tion group for the effective potential, while its value can be 
consistently determined order by order, making it a function 
of the coupling g 2 . It reflects the vacuum energy divergence 
°< j 2 . We refer to [3] for the seminal paper plus toolbox. For 
the current work, we shall not need the explicit effective po- 
tential. As an important asset of Ad quantum field theory is 
its multiplicative renormalizability — to have a clean-cut way 
to deal with UV infinities — , this property needs to be estab- 
lished for eq. (2). Using a more general set of sources, leading 
to an extensive list of Ward identities to be obeyed by the 1PI 
generating functional including the LCOs, this can be proven 
to all orders of perturbation theory [1], from where [24] it 
also becomes evident that there is no pure vacuum term in J. 
Having assured ourselves of the controllable UV nature of our 
effective action, we still need a workable version of it. We em- 
ploy a classical trick: we introduce 2 scalar fields via a double 
Hubbard-Stratonovich (HS) unity into the partition function, 

\ = 9if [ d0t / E ] e -/dVa-O'+O 2 ) 2 e - 2 L / dVE-A/+O l) 2 ) (3) 
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which leads to the following equivalent action, instead of (2), 



s,+s 



F 2 +yIpy-Xd 2 ^-t,d 2 X 



-r\d 2 Q + Qd\ 



2C C 



0? 



+ 



of 



+ -Oi 



2Z, 2A 2A A 



^ d 4 x(-a/-E./ + A/ 2 ) 



(4) 



Here A is a mass dimension 1 parameter, the introduction of 
which is necessary to end up with the appropriate mass di- 
mensions throughout, given that dim[o] = 2, dim[E] = 5/2. 
Notice that A will not enter any physical result at the end if 
we were to compute exactly, as indeed the underlying trans- 
formation (3) constitutes nothing more than a unity. Acting 



with 



5177} 



=/=o 



on both versions of the partition function, 



i.e. before and after the HS transformation, simply provides 
with the correspondences (E) = — (Oi) and (a) = — ((h). 

Let us now collect a few important observations, (i) In 
a loop expansion approach, the parameter A will unavoid- 
ably enter any calculated quantity. However, under the as- 
sumption that passing to higher and higher order one gets 
closer and closer to the exact result, which encompasses A- 
independence, we can fix A in a case-by-case scenario by re- 
lying on the principle of minimal sensitivity (PMS) [4]: we 
must look for solutions of = or higher derivatives if the 
latter eq. has no zeros, (ii) The effective potential itself could 
be examined using the background field method [5]; then we 
do no longer need the sources and can set them to zero, i.e. we 
can drop the Sjj part in the action (4). (iii) If the dynamics 
of the QCD action would decide that (E) = (a) = 0, then we 
are again dealing with nothing else than QCD without D%SB, 
as the trivial and thus irrelevant unities can be integrated out. 
If, on the contrary, (E) ^ by means of dimensional trans- 
mutation, we find ourselves in a vacuum where chiral invari- 
ance is dynamically broken. We draw attention to the crucial 
invariance of the action itself, h$Sf = 0, since we naturally 
have 85E = — 85O1. The role of (a) is to furnish the dynami- 
cal mass for the auxiliary fermion fields, the "m 2 " of [1]. (iv) 
One might wonder about the connection with the standard chi- 
ral condensate, (xjA|/)? This can be straightforwardly inferred 
to be (\j7\|/) oc (E), for example by adding a mass term p to \jn|/ 
in (4) and deriving the vacuum energy w.r.t. p at the end whilst 
setting p = 0. This illustrates clearly the 1-1 correspondence 
between (E), (\j7\|/) and chiral symmetry breaking, (v) With a 
bare quark mass p, the tree level quark propagator yields 



(W>, 



-M{p 2 ) 



M 2 { P 2 



M(p 2 ) = 



< 2 + (a) 



(5) 



One should observe the momentum dependent mass function 
9rf(p 2 ) as a result of the chiral symmetry breaking. Precisely 
this tree level functional form has been applied in fits to lattice 
Landau gauge quark propagators in [7, 8]. These observations 
corroborate the relevance of our model, since we end up with 
a quark propagator in a nontrivial vacuum, which functional 
form is consistent with the nonperturbative lattice counterpart. 



Moreover, we did not make any sacrifices w.r.t. the renor- 
malizability and in the absence of condensation, our action 
is exactly equivalent to QCD. For further usage, we will take 
advantage of these lattice studies and employ the fitting pa- 
rameters to fix a value of our condensates. This will provide 
us with a tree level quark propagator, with the global form 
factor Z = 1 . Including loop corrections on top of the non- 
perturbative vacuum will lead to a nontrivial Z(p 2 ), as also 
seen in e.g. FIG. 5 of [8], which however only deviates mildly 
from 1 over a large range of momenta. The tree approxima- 
tion Z = 1 thus appears to be a valid one. We shall hence not 
attempt here to self-consistently derive and minimize the ef- 
fective action to compute (E) and (a) , which is a technically 
challenging task beyond the scope of this work [3]. (vi) As 
a final important remark, let us scrutinize if we can introduce 
a pion field. We take (E) ^ and consider the already intro- 
duced field Jt [we assume the chiral limit here, p = 0]. Using 
Sf, the chiral current is easily derived to be = vj/y^/y 5 !)/. We 
can then follow more or less the standard derivation, presented 
in e.g. [6]: the correlation function Q^(x — y) = (ft(x)n{y)) 
is subject to d^§ M (x — y) = 8(x — y)(E), as can be shown 
by either using the path integral or a current algebra argu- 
ment as in [6]. Fourier transforming and using the Euclidean 
rotational invariance learns that = 4f (E). To close the 
argument, we consider the 5-matrix element of the current 
destroying a pion state, (£(x)n(p)) « p^ ,px which is re- 
lated to the (pion amputated) propagator when the pion is 
put on-shell. Assuming the pion has a mass m 2 , we would 
get (f 5 (x)%(p)) - \im p2 ^_ m 2(p 2 +m 2 % )g fl (p)e- i P x with pion 

propagator (nn)^_ m i « ^ 
ately provides us with lim /)2 _ 



2 . Recombination immedi- 

K 

nl(P 2 + m l)Qn{P) <* Pn> and 

since we have already shown that Q u = pf (E) , we must re- 
quire m 2 — 0, i.e. the field Jt does describe a massless parti- 
cle. It is instructive to notice that, using the tree level action 
stemming from (4) in the condensed phase, we can rewrite 
71 ^ ¥ _^ / ^ g ^ 7 5 ¥ by means of the equations of motions of the 
auxiliary fermions. We recognize here a nonlocal version of 
the usual pseudoscalar pion field. This is not a surprise, since 
the action (4) itself becomes a nonlocal quark action upon in- 
tegrating out the extra fermion fields in the condensed phase. 

As the proposed framework displays the desired chiral 
properties, we may address further its meson spectrum. As 
a representative example, let us construct and solve a gap 
equation for the (charged) meson mass under suitable 
simplifying approximations. Using the data of [8], in case 
their FIG. 5, we consider degenerate up (m) and down (d) 
quarks with current mass p = 0.014 GeV. The dynamical 
quark mass can be fitted excellently with M(p 2 ) = p 2 M +m2 +P 

with M 3 = 0.1960(84) GeV 3 , m 2 = 0.639(46) GeV 2 with 
X 2 /d.o.f. = 1.18, see FIG. 1. We need = uy M d, p+ = dy^u 
with correlator 



(P„ Pv >* 



>flV 



(6) 



keeping in mind that in the degenerate case, it is transverse 
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thanks to the EOMs and therefore guaranteed to describe the 
3 polarizations of a massive spin 1 particle. We are ultimately 
interested in obtaining a pole in the above channel, corre- 
sponding to a bound state. We rely on the technology set out 
in [9], adapted to the QCD case. Specifically, we shall first 
be concerned with the one-loop contribution to the correlation 
function, 



{2nf 



)4 -r + (*,9)/[(*-9) 2 {{k-qf+m 2 f 

+ (M 3 + p ((* - qf + m 2 )) 2 ] [q 2 (q 2 + m 2 f + (M 3 + 

H(q 2 +m 2 )) 2 ] , f- + (k,q) = (q-(k-q))((k-q) 2 + m 2 ) 2 

X (q 2 + m 2 ) 2 +2(M 3 +/j((k-q) 2 + m 2 )) {{k~q) 2 +m 2 ) 
x(M 3 +p(q 2 + m 2 ))(q 2 + m 2 ) (7) 

We can decompose the appearing denominator into R/(q 2 + 
<a)+R+ / (q 2 + CO, + z'8)+/?_ / (q 2 + C0 r - /8) by using the poles 
3'o = -co, y± = — co r ± i6 = — C0± of the cubic equation 
y(y + m 2 ) 2 + (M 3 +p(y + m 2 )) 2 = 0. In appropriate GeV 
units: R ~ 2.467, co ~ 0.849, R± w -1.234 ± /15.121, co± = 
0.214 ± /0. 052. With the previous numbers filled in, we en- 
counter a real pole and a pair of cc poles, which correspond to 
the so-called /-particles of the seminal works [10] where it has 
been discussed how a pair of such poles can be combined to 
give a physical, i.e. consistent with the Kallen-Lehmann (KL) 
representation, contribution to the spectral form of a bound 
state propagator [25]. For the remainder of this work, we shall 
thence only be concerned by this physical part of the correla- 
tor, under the assumption that any unphysical piece, originat- 
ing from combining poles in pairs that are not cc, will eventu- 
ally cancel out when a full-fledged analysis and methodology 
to deal with the /-particles would become feasible. Explicitly, 
we can extract the spectral representation using the complex 
mass cut-rule technology of [11]. A rather tedious computa- 
tion leads to 



(?a?a)k= / d ^ 



(X) 



o x + k 2 
R' 



unphysical , p + (x) > 



p~ + (x) = 9(x-Xi)^V 1 /4-co/x (x/2-co)(co + m 2 



2\4 



-2(M 3 +p((o + M 2 )) (co + m 
v/(x-2co r ) 2 -4(co 2 + e 2 ) 



-9(x-x 2 )^ 



(x/2-co,.)((co r + m 2 ) 2 + 9 2 ) 2 



+2((M 3 + i u(co r + m 2 )) 2 + i u 2 9 2 ) ((co r + ; 



,2\2 



-e 2 



(8) 



with thresholds X] = 4co, x 2 = 2co r + 2 v /co 2 T9 2 . Eq. (8) is 
the spectral representation of the single bubble approximation 
to the p-correlator. In order to find a bound state, we need to 
take into account the QCD interaction structure. The quarks 
do not interact directly, the force being mediated by the gluon. 
Unfortunately, taking the full gluon interaction into account 
appears to be a hopelessly complicated task. To simplify the 
analysis to get a tractable problem in this first attempt, we 




FIG. 1: Lattice quark mass function [8] with its fit (M(p 2 ). 
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FIG. 2: Bubble diagrams for the meson correlator. Full (open) circles 
represent four-fermion vertices (meson operator insertions); solid 
lines are nonperturbative quark propagators. 



consider the observation made in e.g. [12] that a gluon contact 
point interaction can give qualitatively good results. Specifi- 
cally, it is by now well accepted that the Landau gauge gluon 
propagator (D(p 2 ) becomes dynamically massive-like in the 
infrared. We make the assumption that the gluon is massive- 
like in the region p < 1-2 GeV, which is anyhow the region 
where the relevant QCD physics is supposed to happen and 
that it can be approximated by a constant in momentum space, 
2)(/7 2 ) = 4f, or by (D(x ~ y) « 8(x — y) in position space. In- 
tegrating out at lowest order such gluon leads, via the Dirac 
delta behaviour, to a NJL-like (contact) interaction between 
the quarks, more precisely one finds after some Fierz rear- 
ranging, see the Appendix of [13], the effective interaction 
\G{q+q) 2 + lower in \/N, where ♦ e {1,/y 5 , -^y A „ -^7^75} 



and G = 



A 2 N 



To allow for additional simplification with- 



out throwing away too much crucial dynamics, we may fur- 
thermore consider only the leading order in 1 /N. For the 
eventual effective coupling G, we could borrow the various 
NJL estimates from [14], where G ~ 6-10 GeV -2 . We can 
however also estimate more directly an appropriate value. 
In the Landau gauge, the gluon and ghost propagator form 
factors, Z g i(p 2 ,p 2 ) and Z g f,(p 2 ,~p 2 ), can be used to define 
a renormalization scale p independent strong coupling con- 
stant, 4ng 2 { P 2 ) = a( P 2 ) = a(p 2 )Z gl (p 2 ,il 2 )Z 2 h (p 2 ,p 2 ), see 
e.g. [15]. Using the most recent available lattice data on this 
matter [16], albeit for = 2 without fermions, which rely on a 
MOM scheme (Z g h = Z g i = 1 at p 2 = p 2 , with p = 2.2 GeV), 
we can estimate from their FIG. 5 that a(p) ~ 0.5. For the 
"constant" MOM scheme gluon propagator, we may sim- 
ply set A 2 = p 2 , thereby overestimating the UV and under- 
estimating the IR, to roughly approximate G ~ 8.5 GeV -2 , 
quite nicely in the NJL ballpark. For this paper, we will use 
G = 5,7.5,10 GeV -2 as exemplary values. In the large-Af 
approximation, we can then consistently consider the sum of 
bubble diagrams [26] in our quark model, see FIG. 2. 

The four-fermion coupling includes a priori all interaction 
channels ♦ 6 {MY) TTfY/J! 'J2 ^ 5 }' wn ^ e ^ e initial and final 
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insertions carry the vectorial character of the p meson. In this 
case, as well-known, only the vector four-fermion coupling 
Y„ contributes [27] and the resummation of the bubble chain 
reduces to a geometric series involving the one-loop result, 
(8). Using the already derived spectral form, we then get 



■K.- + (k 2 ) = 



?- + {k 2 ) 
l + l!F-+(k 2 ) 



■(x)dx 
+ k 2 ' 



An important observation is that the branch cut structure of 
H{_ v {k 2 ) is completely determined by that of f v {k 2 ), in 
particular %_ v (k 2 ) has a physical KL form if T~ + {k 2 ) does. 
The remaining task at hand is to identify if v {k 2 ) allows 
for poles in the physical region max(— Xi,— X2) < k 2 < 0. 
Therefore, we need to solve the gap equation f v {k 2 ) = 
— 2/G. As common for Ad quantum field theories, the quantity 
f v (k 2 ) is divergent due to violent UV behaviour, explicitly 
visible from the integral representation (8). We assume the use 
of dimensional regularization and Landau gauge renormaliza- 
tion factors to kill off sub-loop divergences, whereas the resid- 
ual infinities in the composite operator Green function can be 
taken care off by additive subtractions in the BPHZ approach. 
At the level of the KL integral (8), n > subtractions at scale 
1 > max(— Xj ,— X2) corresponds to 



^ ub + (k 2 ,<T) = (>r-k 2 y 



p-+(x) 



(x + T)«(x + /t 2 )' W 



In the current case at least 2 subtractions are required. Notice 
that if no subtractions were to be necessary, then f- + {k 2 ) 
is a strictly decreasing function thanks to p~ + (x) > 0, so at 
most one solution is possible. This property is not necessarily 
maintained at the subtracted level, except for n = 1 . Conse- 
quently, sometimes spurious extra solutions appear, caused by 
the enforced functional behaviour after subtraction. For this 
reason, we take n = 3 in which case only a single solution is 
found. Next to the pole of the propagator, also the associated 
residue carries physical information. With the conventions of 
[18, 19], we define the decay constant / p ± via /p±m p ±e ;U = 
(0|7Jy ;U af|p + ), with the polarization tensor of the p + meson, 
normalized as £ M • £u = 3. Using the matrix element represen- 
tation of the KL spectral density, it becomes clear that 3f 2 ± mp 
is nothing else than the residue of H{_ v {k 2 ) at its pole —m 2 ± . 
In FIG. 3, we have displayed both m p (T) and / p ±(T). To 
get a grip on a reasonable value for T, we rely again on PMS 
as observable quantities should not depend on a chosen sub- 
traction scale. Since we have two (related) physical quantities 
at hand, we minimized S(T) = |ra p± (T)| + |/ p ±(T)| where 
the prime means derivation w.r.t T, and we used the rescaled 
mass and decay constant, ra p ±(T) = m p ±('T)/m p ±(— X2), 

fp ± ( r T) — fp ± ( r ^)/fp ± (~^2) to attribute to both quantities 
an "equal start". We choose the smallest possible scale, 



viz. 1 = — X2, as reference. FIG. 3 shows that 8(T) develops a 
clear minimum for all values of the effective coupling. Specif- 
ically, we find % G=5 ' 7 5 ' 10 w -0.38, -0.43,-0.46 GeV lead- 



G=5, 7.5, 10 



i 0.84,0.83,0.83 GeV and f^T 5 ' 7 ' 5 ' 10 : 



ing to m p± 

0.13,0.10,0.09 GeV. We notice that the mass estimate is 
pretty stable, while the decay constant seems to be more sen- 
sitive to the value of the coupling. We found it useless to 
attempt presenting a detailed error analysis given the vari- 
ous deliberately rude approximations, but we find it at least 
reassuring that we do find a result not too far off the ex- 
perimental p meson mass, 775.49 ±0.34 MeV [17], despite 
the approximations made, next to having used heavier bare 
quarks. For the decay constant, we can quote experimental 
and lattice estimates, f p f « 0.208 GeV [20], f l p f< « 0.25 GeV 
[18]. Our estimates, rude as they may be, are somewhat 
lower, though it must be remarked that the heavier the particle 
mass, the smaller the decay will get for a fixed value of the 
residue. If we, for example, fix T to its value at G — 5 GeV 2 
where m p ± equals its experimental value, we would obtain 
f p ± rj 0. 16 GeV, whose deviation from the experimental value 
is of the same size as its lattice counterpart. 

Summarizing, we presented a calculable quark model with 
the right properties to describe the mesonic QCD sector. In fu- 
ture work, it would, for example, be interesting to (i) compute 
the effective potential to get self-consistent estimates for the 
condensates, (ii) consider other meson states and their prop- 
erties, (iii) derive the Gell-Mann-Oakes-Renner like relations 
for multiflavour versions with bare quark masses, etc. Given 
that our model incorporates an effective way of confining 
quarks simultaneously with breaking the chiral symmetry, 
thereby generating propagators that are consistent with lattice 
QCD, it might provide an interesting setup to investigate non- 
trivial finite temperature dynamics [21], given the intertwining 
of (de)confinement and chiral symmetry breaking/restoration. 




G = 5 GeV- 2 



C = 7-5 GeV - 



-Q.5 0.0 0.5 




FIG. 3: m„± (full line) and f p ± (dashed line) in terms of 1 (left) and 
8(T) (right) (see main text for definitions). 
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